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We investigate the effect of spin-orbit coupling on the band structure of graphene-based two- 
dimensional Dirac fermion gases in the quantum Hall regime. Taking monolayer graphene as our 
first candidate, we show that a quantum phase transition between two distinct topological states - 
the quantum Hall and the quantum spin Hall phases - can be driven by simply tuning the Fermi 
level with a gate voltage. This transition is characterized by the existence of a chiral spin-polarized 
edge state propagating along the interface separating the two topological phases. We then apply 
our analysis to the more difficult case of bilayer graphene. Unlike in monolayer graphene, spin- 
orbit coupling by itself has indeed been predicted to be unsuccessful in driving bilayer graphene 
into a topological phase, due to the existence of an even number of pairs of spin-polarized edge 
states. While we show that this remains the case in the quantum Hall regime, we point out that by 
additionally breaking the layer inversion symmetry, a non-trivial quantum spin Hall phase can re- 
emerge in bilayer graphene at low energy. We consider two different symmetry-breaking mechanisms: 
inducing spin-orbit coupling only in the upper layer, and applying a perpendicular electric field. In 
both cases, the presence at low energy of an odd number of pairs of edge states can be driven by an 
exchange field. The related situation in trilayer graphene is also discussed. 

PACS numbers: 72.80.Vp, 73.43.Nq, 73.20.At, 73.21.Ac 



I. INTRODUCTION 

Topological insulators are bulk insulators which pos- 
sess robust conducting surface states [1, 2], Paradigmatic 
two-dimensional examples of this class are the quantum 
Hall (QH) and the quantum spin Hall (QSH) phases, 
which are characterized by respectively chiral and he- 
lical one-dimensional edge states. While the former can 
be generated by simply applying a strong perpendicular 
magnetic field and has been rather extensively studied 
since the 1980s, the latter requires the presence of spin- 
orbit coupling and has received very little experimental 
evidence. Indeed, despite the wide interest shown in the 
literature for the QSH phase (and for topological phases 
in general) since the seminal works by Kane and Mele 
[3, 4], experimental traces of this phase have remained 
scarce, with the exception of the remarkable works on 
HgTe quantum wells [5-7] (see also the experiment in- 
volving In As [8]). Recent studies [9-11] have revived 
the possibility of generating a QSH in graphene [12, 13], 
by showing that low concentrations of suitably chosen 
adatoms, randomly deposited on graphene, could open a 
large non-trivial gap in graphene's otherwise semimctallic 
band structure, and yield transport properties showing 
no trace of the spatially inhomogenous spin-orbit cou- 
pling. The perspective of successfully turning graphene 
into a QSH insulator is promising, as it would consider- 
ably enhance the experimental feasibility of engineering 
samples of the latter, which are so far limited to the pre- 
viously mentioned and experimentally challenging HgTe 
heterostructures. 

Although it does not enjoy the conceptual simplicity of 
the monolayer, bilayer graphene is an interesting system 
in its own right. It is a gapless semimetal, characterized 
by massive chiral excitations carrying a topological Berry 



phase 2ir [14], with a very rich list of many-body insta- 
bilities predicted at low density (see [15] for correspond- 
ing references). One of its most remarkable properties, 
in contrast to monolayer graphene, is the possibility to 
open a gap in its band structure by simply applying a per- 
pendicular electric field which breaks the layer inversion 
symmetry [16]. However, in the presence of a perpen- 
dicular magnetic field, the electronic properties of both 
systems become qualitatively very similar. In particular, 
their energy spectrum is characterized by a particle-hole 
symmetry and by the existence of levels sitting exactly 
at zero energy which are at the origin of the anomalous 
quantization of the Hall conductance in these systems as 
compared to other two-dimensional electron gases [17— 
19]. This is the hallmark property of what we will refer 
to in this article as two-dimensional Dirac fermion gases 
(2DDFGs). 

The purpose of this article is to discuss what happens 
to the band structure of a 2DDFG when the effects of 
both magnetic field and spin-orbit coupling are taken into 
account simultaneously. Taking graphene as our first ex- 
ample, we shall review in section II the results already 
published elsewhere [20], according to which a topologi- 
cal phase transition takes place at low energy and can be 
tuned by simply varying the chemical potential. Then, 
in the following section, we shall investigate the related 
situation in bilayer graphene and show that the results 
obtained for the monolayer do not extend to it. This can 
be traced back to the fact that the topological invariant 
characterizing the QSH phase is non-trivial only if there 
are an odd number of pairs of edge states, which trans- 
lates in multi-layer graphene into the condition of having 
an odd number of layers. Nevertheless, we will show in 
section IV that by adding additional ingredients to our 
model, namely by breaking the layer inversion symme- 
try, a non-trivial QSH phase can be generated in bilayer 
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graphene along with a corresponding topological phase 
transition. We stress that (almost) all of our results can 
be understood by simply looking at the band structures 
of the systems we study. Finally, section V discusses the 
possible extension of our approach to other systems, and 
we conclude in section VI. 



II. GRAPHENE 




Recent investigations of the interplay between QH and 
QSH phases in some specific examples of 2DDFGs [20- 
23] have led to surprising results. In these works, it was 
shown that the QSH phase can survive the presence of 
a perpendicular magnetic field and that the Z2 topolog- 
ical invariant [4] remains non-trivial for energies below 
the spin-orbit induced gap, despite the breaking of time- 
reversal symmetry [59] . As we will exemplify in the case 
of graphene, the origin of this intriguing result actually 
stems from the existence of zero-energy Landau levels: 
as soon as the spin-degeneracy of these levels is lifted, 
spin-polarized edge states characteristic of a QSH phase 
emerge [60] . 



A. Model 

Let us start by introducing the model from which our 
results shall be derived. In the vicinity of the zero-energy 
points in the Brillouin zone, low-energy excitations can 
be described by a Dirac Hamiltonian: 



Hg = v F (Tp x a x + Py(Ty) , 



(1) 



with a x , <Ty the usual set of Pauli matrices acting in the 
two-dimensional space of the inequivalent sublattices A 
and B (see Fig. 1), and with vp = \/3toa/(2h) the Fermi 
velocity, expressed as a function of the microscopic lattice 
parameters to (nearest-neighbor hopping amplitude) and 
a (lattice constant) which we choose in the following as 
our working units of energy and length, respectively, r = 
±1 accounts for the two possible valleys from which the 
low-energy excitations can arise. In the vicinity of these 
points, the energy-dispersion relation reads e = ±«j?|p|. 

The presence of a perpendicular magnetic field can be 
straightforwardly included by making use of the Peierls 
substitution p > II = p + eA, which accounts for the 
presence of the magnetic vector potential A such that 
Vx A = Bz. The components of the generalized momen- 
tum satisfy the Heisenberg algebra [II x ,II y ] = —i(h/l B ) 2 . 
By expressing these components in terms of the usual 
harmonic oscillator ladder operators [19], 



n, 



-^-(a + ot),n B = -^-(a-ot), (2) 



and using the standard raising and lowering properties 
of these operators on the eigenstates (a\n) = \fn\n — 1) 
and a)\n — 1) = y/n\n)), the energy spectrum can then 



Figure 1: Sketch of the graphene hexagonal lattice, character- 
ized by nearest-neighbor hopping to and next-nearest-neigbor 
Kane-Mele spin-orbit coupling A so . The unit cell of the lattice 
contains two inequivalent sites, A and B. 



straightforwardly be shown to turn into the well-known 
Landau levels, 



(3) 



with Ab = \f2twp jls , and Is — \Jh/{eB) the magnetic 
length. As already mentioned before, the main distinc- 
tive feature of the Landau level spectrum of a 2DDFG as 
compared to that of a standard two-dimensional electron 
gas is the existence of a zero-energy level at n = orig- 
inating from the pseudo-relativistic nature of the charge 
carriers. Also note that all levels enjoy a 4-fold degener- 
acy arising from spin and valley indices. 



B. Band structure / edge state correspondence 

The appearance of edge states in this context can be 
best understood by looking at the band structure of a 
graphene ribbon. The latter is a system which is trans- 
lationally invariant in one direction, and confined in the 
other. In order to derive the band structure numerically, 
we formulate the above ingredients in terms of a tight- 
binding model, in which Eq. (1) becomes 



(4) 



Indices label lattice sites, while symbol ( ) refers 

to nearest-neighbor coupling (with hopping amplitude 
to), as is illustrated in Fig. I. The Peierls phase ^ij — 
(e/fi) J" A • dr takes into account the contribution from 
the magnetic flux threading the lattice. Numerical cal- 
culations are performed using kwant, the new quantum 
transport software package developed by A. Akhmerov, 
C. Groth, X. Waintal, and M. Wimmer. In the process, 
we choose to work with armchair boundary conditions, 
but our results are qualitatively unaffected by this choice. 

The band structure associated with the peculiar spec- 
trum of Eq. (3) in a ribbon geometry (translationally 
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Figure 2: (Color online): Energy spectrum of a monolayer 
graphene armchair ribbon in the QH regime (C / 0). Black 
circles and red triangles respectively stand for spin up and 
spin down bands, which are here indistinguishable due to spin 
degeneracy. The horizontal dashed line represents an arbitrar- 
ily chosen Fermi level. Its intersection with the lowest band 
of the system indicates the existence of a (spin-degenerate) 
edge state which propagates in opposite directions on oppo- 
site sides. The ribbon width is W = 38 and the magnetic 
length is Ib ~ 4. 



invariant in the x-direction, confined in the y-direction 
with \y\ < W/2) is shown in Fig. 2. Due to the nature 
of the classical dynamics, the transverse coordinate of 
the cyclotronic center of motion y c can be identified with 
the conserved longitudinal momentum k x via the formula 
y c = —k x l 2 B . Observing the band structure in terms of 
this real space coordinate, one can see in Fig. 2 that, far 
from the edges of the ribbon, the band structure consists 
of flat bands which are none other than the Landau levels 
of Eq. (3): electrons in the bulk are classically localized 
by the magnetic field along closed cyclotronic orbits. On 
the other hand, electrons in the vicinity of the edges can 
scatter along them and propagate following skipping or- 
bits, which translates in the band structure into bulk 
Landau levels acquiring a finite dispersion as they ap- 
proach the edges of the ribbon: 



,(n) _ 



1 de n 
h dk x 



(5) 



Because this dispersion is monotonous on a given edge 
(see Fig. 2), the edge states cannot be backscattered un- 
less they are coupled to the states living on the opposite 
edge, a process the likelihood of which decays exponen- 
tially with the width of the system. This property of 
the edge states is generally referred to as chirality and 
is the reason why these states can carry current without 
dissipation: this leads to the celebrated QH effect [24], 
characterized by a quantized conductance G = G(e 2 /h). 
More formally, the edge states enjoy a topological pro- 



tection encoded in the Chern number C which is a Z 
topological invariant characterizing the number of filled 
bands in the QH regime [25] . It is a topological quantity, 
in the sense that smooth deformations of the Hamilto- 
nian (deformations which do not close the gap) cannot 
change its value, and shall be defined in the next subsec- 
tion. We thus see that, for most purposes, the physics 
of topological phases such as the QH phase can be very 
simply extracted from the corresponding band structure. 

We now consider the situation where, in addition to 
the perpendicular magnetic field, the effect of spin-orbit 
coupling as introduced by Kane and Mele [3] is accounted 
for in the Hamiltonian as 

H so = rsA so a z , (6) 

which is characterized by the energy scale A so . r = ±1 
and s = ±1 account for valley and spin degrees of free- 
dom. In terms of a tight-binding model, Eq. (6) can be 
implemented as [3] 



(7) 



where indices once more label lattice sites, while 

(a, (3) label spin indices, symbol (( )) refers to next- 
nearest-neighbor coupling (with SO-induced hopping am- 
plitude A so = A so / (3\/3) [3]), and i/y = ±1 depending on 
whether sites are coupled clockwise or counter-clockwise 
(see Fig. 1). Note that in order for the system to remain 
gauge invariant, Peierls substitution has to be done on 
all hopping matrix elements: nearest-neighbor and (SO) 
second nearest-neighbor. The presence of spin-orbit cou- 
pling modifies the Landau level spectrum according to 
the expression 



±a/A 



f A s 2 , for n ? 
-sA so , for n = 



(8) 



The latter is characterized by the n = level being lifted 
from zero energy into spin-polarized branches: E = +A so 
features only spin-down states, while E = — A so features 
only spin- up states [22] (see Fig. 3). While other Lan- 
dau levels retain their associated chiral edge states, ir- 
respective of the spin polarization, the lowest Landau 
level now features counter-propagating edge states for 
\E\ < A so , with a spin-dependent direction of propa- 
gation (see Fig. 3): one has v x '^ ■ v x '^ < on a given 
edge (allowing for an additional spin dependence in the 
definition of Eq. (5)). This is illustrated in the real-space 
sketch above the band structure in Fig. 3. It is the signa- 
ture of a QSH phase, as can be certified by computing the 
associated Z2 topological invariant introduced by Kane 
and Mele [4], which we do next. 



C. Topological order 

Let us start by recalling the standard topological num- 
ber characterization of Landau levels when A so = 0. 
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Figure 3: (Color online): Same as in Fig. 2, but with an ad- 
ditional spin-orbit coupling term A so = 0.02. The latter lifts 
the spin-degeneracy of the zero-energy Landau level, yielding 
a QSH phase (y = 1) with a single pair of counter-propagating 
spin-polarized edge states. The sketch above the band struc- 
ture depicts the edge states in real-space (thick black vertical 
lines on the sides represent the edges of the ribbon). 



Each Landau level n and its associated eigenfunctions 
over the first Brillouin zone are characterized by a topo- 
logical invariant, the so-called Chern number [25]. This 

topological number takes a value Cr"i = +1 for each Lan- 
dau level, independently of the Landau n, valley r or 
spin s indices. For each value of the Fermi energy, we 
can characterize the corresponding phase by a topologi- 
cal number 



e T , s , with e T , s (E F ) = J2 &( r 



(9) 



obtained by summing the Chern numbers of all filled en- 
ergy bands [25]. 

In 2DDFGs, the Chern number is a priori ill-defined 
because of the existence of an infinite number of filled 
energy bands of negative energy. Through the use of non- 
commutative Berry's connection, it was however shown 
[26] that the Chern number takes a value Q T:S {E = _ ) = 
— 1/2 per degree of freedom for energies immediately be- 
low the Dirac point. In the case of graphene, defining the 
spin Chern number as 



£e T , 



(10) 



one obtains G S (E = ) = — 1 per spin species (since 
there are two valleys). With this prescription, the band 



structure of graphene in the QH regime can be easily 
described by computing the value of the Chern number 
as a function of the Fermi energy. This yields 



e t 
e t 



-1 , for — A.b < Ep < 
2n + 1 , for E F > 



(11) 



where n is the index of the highest filled Landau level. 
Notice that Chern numbers of each spin species are equal 
(since the spectrum is spin-degenerate) and that 6 = 
G-f + Cj, is always non-zero, as expected for a QH phase. 
The total Chern number increases step-wise by multiples 
of 4, due to spin and valley degeneracy of the Landau 
levels. 

Restoring a finite value of A so , the topologically non- 
trivial nature of the phase for E < A so can be checked 
by computing the corresponding value of the Z2 topo- 
logical invariant. In the presence of spin rotational sym- 
metry (conservation of S z ), this invariant can be simply 
expressed as the difference of the Chern numbers for each 
spin species [27]: 



*, = i(e t -e+) (mod 2) 



(12) 



Note, however, that the existence of this invariant is nat- 
urally independent of the existence of this symmetry. Fol- 
lowing the calculations performed in [20] , one easily finds 



e t = -64. = 
e t = e l = 2r 



-1 , for \E F \ < A sc 
+ 1 , for E F > A s 



(13) 



where n is once more the index of the highest filled Lan- 
dau level. This time, one is faced with a QH phase for 
\Ep\ > A so , characterized by the same Chern number 
as in the previous case, while for \Ep\ < A so , the total 
Chern number vanishes 6^ + 6^ = 0, indicating that this 
region is no longer in the QH phase. However, as the Z2 
invariant v = 1 does not vanish, the phase in this region 
is a QSH phase. 

The resulting band structure is thus particularly inter- 
esting: it consists of a QSH phase at energies \E\ < A so 
and a QH phase at other energies. One can therefore 
observe a topological phase transition in this system by 
simply tuning the Fermi level across the spin-orbit gap. 
The central manifestation of this phase transition is the 
existence of a spin-polarized state localized at the in- 
terface between both phases. This state could be most 
clearly observed experimentally by making use of an ad- 
ditional electric gate to independently tune the Fermi 
level in two different parts of the system, thereby realiz- 
ing a topological heterojunction. We refer the interested 
reader to our previous paper [20] for a detailed discussion 
of these matters. 



III. BILAYER GRAPHENE 

We now switch to the slightly more involved case of 
(Bernal-stacked) bilayer graphene and start by address- 
ing the possibility of inducing a QSH phase in bilayer 
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Figure 4: Side view of Bernal-stacked bilayer graphene, which 
is characterized by intra-layer nearest-neighbor hopping to 
and inter-layer hopping t\. 



graphene. This is not a trivial endeavour, as the naive 
extension of the Kane-Mele model to bilayer graphene 
yields a weak Z2 topological phase, characterized by an 
even number (rather than an odd number, as in mono- 
layer graphene) of pairs of spin-polarized edge states 
[28, 29]. This doubling of the number of edge states ba- 
sically arises because a bilayer has twice as many layers 
as a monolayer. For the same reason, a graphene trilayer 
will have an odd number of pairs of edge states and there- 
fore feature a non-trivial QSH phase. Breaking the layer 
symmetry by considering the case where spin-orbit cou- 
pling is present in only one of the layers was shown not 
to be any more effective [28] , the system then remaining 
semimetallic. 

Here we follow a different approach, inspired by the 
model presented in the previous section. This seems like a 
natural idea, as bilayer graphene is also known to feature 
zero-energy Landau levels [14, 30]. In this section, we will 
show that the presence of both spin-orbit coupling and a 
perpendicular magnetic field in the bilayer yields a band 
structure very similar to that of monolayer graphene, but 
with results no different from that of Prada et al. [28]: 
the obtained QSH phase is topologically trivial due to 
the existence of an even number of pairs of spin-polarized 
edge states. On the other hand, we will show in the next 
section that if spin-orbit coupling is present in only one 
of the two layers or if a perpendicular electric field is 
applied, then the breaking of layer inversion symmetry 
opens the door for a non-trivial QSH phase to arise at 
low-energy. 



A. Model 

The Hamiltonian for bilayer graphene can be expressed 
using two sets of Pauli matrices {a, rj} which respectively 
refer to sublattice (A, B) and layer (1,2) spaces. We con- 
sider the usual Bernal stacking, inherited from graphite 
(see Fig. 4), in which A atoms in the upper layer (2) lie 
above B atoms from the lower layer (1). Starting from 
the basis (A2, B2, Ax, Bi) T , the Hamiltonian reads [30]: 

Hbg = v F {rp x a x + p y a y )r/o + ^-(a x r) x - a y r) y ) . (14) 



The first term describes the usual low-energy Dirac struc- 
ture of monolayer graphene. The second term takes into 
account the coupling between both layers, characterized 
by an energy scale t\ ~ 0.15. Corrections to Eq. (14) 
such as trigonal warping are small effects, typically only 
relevant below the meV range [15], and will therefore be 
neglected. 

The spectrum associated with Eq. (14) is particle- 
"hole" symmetric, with high-energy bands at e hlgh = ±ii 
and low-energy bands touching at two Dirac points char- 
acterized by a topological Berry phase 27r. In the vicinity 
of this point, the energy-dispersion relation is quadratic, 
e low = ±p 2 /(2m*), with m* = t 1 /(2v 2 F ) the effective 
mass of the gapless excitations. 

The presence of a perpendicular magnetic field can be 
straightforwardly included by making use of the Peierls 
substitution as before, and the energy spectrum can then 
be shown to turn into the well-known Landau levels [30] , 

e n = ±hui c ^n(n - 1) (15) 

with lo c = eB/m* the characteristic cyclotronic fre- 
quency. The latter can be related to the monolayer 
graphene energy scale by the simple relation fku c — 
Ag/ii. Notice how the spectrum in Eq. (15) fea- 
tures twice as many zero-energy levels as in monolayer 
graphene, since the n = 1 level also vanishes. Actually, 
one can prove on general grounds that chirally stacked 
iV-layer graphene should feature a 4iV-fold degenerate 
zero-energy Landau level [31, 32]. One should also have 
in mind that the accuracy of expression (15) for n^0,l 
is only correct in the limit "C t\. 



B. Quantum Hall regime 

To compute the associated band structure numerically, 
we once more formulate the above ingredients in terms 
of a tight-binding model, in which Eq. (14) becomes 

X BG = -t J2e^4c J +t 1 ]T c\c J} (16) 

(ieA 2 ,j£B 1 ) 

using the same notations as in the previous section. No 
Peierls phase appears in the second term, as A 2 -Bx bonds 
are oriented along the z-axis. The band structure asso- 
ciated with the Landau level spectrum of Eq. (15) in a 
ribbon-geometry is displayed in the upper panel of Fig. 5. 
As expected, it resembles very closely that of monolayer 
graphene in the QH regime. The main difference between 
the two lies in the existence of twice as many dispersing 
branches arising from the lowest Landau level in bilayer 
graphene, due to the doubling of the zero-energy Lan- 
dau level degeneracy. This translates in the language 
of topological invariants into Chern numbers per spin 
species C s = —2 immediately below zero energy. With 
this prescription, the band structure of bilayer graphene 
in the QH regime can be easily described by computing 
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C. Effect of layer-symmetric spin-orbit coupling 
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Figure 5: (Color online): (Top panel) Energy spectrum of a 
bilayer graphene armchair ribbon in the QH regime. Black cir- 
cles and red triangles respectively stand for spin up and spin 
down bands, which are here indistinguishable due to spin de- 
generacy. Notice that the zero-energy Landau level has twice 
as many bands as its counterpart in monolayer graphene. The 
ribbon width is W = 38 and the magnetic length is is ~ 4. 
(Bottom panel) Same as above, but with an additional spin- 
orbit coupling term A so = 0.02. The latter lifts the spin- 
degeneracy of the zero-energy Landau level, yielding a weak 
QSH phase (y = (mod 2)) with an even number of pairs of 
counter-propagating spin-polarized edge states. 



the value of the Chern number as a function of the Fermi 
energy, yielding: 



-4- 



-2 , for - hu c V2 < E F < 



C T = G l = 2[max(0, n - 1) + 1] , for E F > 



(17) 



where n is the index of the highest filled Landau level. 
Notice that Chern numbers of each spin species are equal 
(since the spectrum is spin-degenerate) and that 6 — 
Ct + ^4- i s always non-zero, as expected for a QH phase. 



We now consider the situation where, in addition to 
the perpendicular magnetic field, the effect of spin-orbit 
coupling as introduced by Kane and Mele for mono- 
layer graphene is accounted for symmetrically in both 
layers [61]. The layer-degenerate Kane-Mele spin-orbit 
coupling term is encoded in the Hamiltonian 



H ao = TsA so a z T] Q 



(18) 



which, in terms of a tight-binding model, can be imple- 
mented as 

?4o = *A so J2 ^ e ^'( c l,a s f c i,0) , (19) 

((id)) 

with similar notations as in the previous section, symbol 
(( )) referring to intra- layer next-nearest-neighbor cou- 
pling. The presence of spin-orbit coupling modifies the 
Landau level spectrum according to the expression 



±y/n(n- l)(fiw c ) 2 + A 2 Q , for n ^ 0, 1 
— sA so , for n = 0, 1 



(20) 



The latter is characterized by n = and n = 1 levels 
lifted from zero energy into spin-polarized branches: E = 
+A so features only spin-down states, while E = — A so 
features only spin-up states (see lower panel of Fig. 5) . 

The topologically trivial nature of the corresponding 
low-energy phase can be checked by computing the value 
of the Z2 topological invariant. As a straightforward gen- 
eralization of the calculations performed in [20], one ob- 
tains the following results: 



e t 



-4- 



= +2 , for \E F \ < A s 
2[max(0,n- 1) + 1] 



for E F > A s 



(21) 

where n is once more the index of the highest filled Lan- 
dau level. This time, one is faced with a QH phase 
for \Ep\ > A so , characterized by the same Chern num- 
ber as in the absence of spin-orbit coupling, while for 
\Ep\ < A so , the total Chern number vanishes C^+C^ = 0, 
indicating that this region is no longer in the QH phase. 
However, as the Z2 invariant v — (mod 2) also vanishes, 
the phase in this region is not a QSH phase either: rather, 
it is a weak QSH phase, in the sense that time-reversal- 
symmetric perturbations can couple the edge states and 
induce backscattering. This was not the case in mono- 
layer graphene, due to the existence in the latter of a sin- 
gle pair of counter-propagating spin-polarized edge states 
at low energy. In the situation discussed in this section, 
we are thus led to conclude that a similar picture as that 
described in Ref. [28] prevails. The way around this in- 
volves breaking the layer inversion symmetry, as we will 
see in the next section. 

Before moving on, however, we would like to pause 
and comment on the fact that the model we considered 
in this section could provide a convenient platform for 
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testing precisely how weak a v — (mod 2) QSH phase 
would be. Indeed, even though theory predicts that pairs 

of edge states should couple through backscattering pro- — 1 — 1 — I — 1 — r~7i — 1 — 1 — 1 — 

cesses, an experimental measure of how strongly edge ". *.". \ Cvp = C-A = ^ / •' 

state transport would be destroyed by such processes is °-l *• ''y ^^fc^^wtw^t^j. «?. . , , 

yet to be done, and one cannot exclude the possibility of • •) _ 

unexpected robustness, similarly to what has recently be- *'"'*•, y^-^=- 2 y^^. = ..Q....... 

gun to be understood in so-called weak three-dimensional •. .Ay""* yr^ <E^^1*|* *^^> • 
topological insulators [33, 34]. Said a little differently, w ' Cf= l C*=-2 ^\ 

it remains unclear how one could distinguish through A * / s-y \ ' 

transport measurements a topological phase from a triv- •* •' = v>^~* 

ial phase which has edge states (such as the one exhibited ■ / ~ (3 f = .j - -2 ~ \ 

in this section). " ' ; /." ' "' """" " ' "" 



IV. BILAYER GRAPHENE WITH BROKEN 
LAYER INVERSION SYMMETRY 
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This section is devoted to the study of two layer in- 
version symmetry breaking mechanisms which enable an 
exchange-induced QSH phase to arise at low energy: (i) 
inducing spin-orbit coupling only in one of the layers, and 
(ii) applying a perpendicular electric field. We provide 
estimations of the magnitude of the QSH gap, and also 
address other possible phases which appear in our set- 
tings: a spin-polarized QH phase and a quantum valley 
Hall phase. 

A. Mechanism I: layer-asymmetric spin-orbit 
coupling 

Using the same basis as in Eq. (14), we now replace 
Eq. (18) by the symmetry-breaking term 

H:r^r S A so a z (^±^^j (22) 

which induces spin-orbit coupling only in the upper layer. 
This situation is physically relevant if one considers the 
possibility of inducing spin-orbit coupling in graphene 
by depositing adatoms on the surface [9-11]. The cor- 
responding tight-binding expression is given by applying 
Eq. (19) only in the upper layer. 

In order to obtain the new Landau level spectrum, us- 
ing the ladder operators introduced in section II, one 
must now solve a quartic equation e„+a n e„+ /3 n e n+7n = 
0, with a non-vanishing linear term /3„ 0: 

a n = -(*? + Al + (2n - 1)A|) 
n = -sA so tj . (23) 

7„ = (n - l)A|(nA| + A S 2 J 

Note that the role of opposite spin polarizations is ex- 
changed when going from the conduction to the valence 
band: e„(— s) = —e n (s). This leaves 4 (out of the 8) 
eigenvalues of the lowest Landau levels to be found. Two 
can easily be identified: cq = — sA so (for n = 0) and 




hco /t, 

c 1 



Figure 6: (Color online): (Top panel) Band structure with 
same parameters as in Fig. 5, except that spin-orbit cou- 
pling is only applied in the upper layer. This time, a 
spin-degenerate zero-energy Landau level survives, while low- 
energy edge states are characterized by an unbalanced spin 
population on a given edge (|Cf| ^ | C4. | ) . (Bottom panel) 
Dependence of the eigenvalues e± on the magnetic field for 
A so = 0.01. Dashed lines correspond to the analytical pre- 
dictions (valid in the perturbative limit hw c <C ti) and thick 
lines to the numerically obtained values. Inset: Dependence 
of e± on the spin-orbit gap for A_b = 0.05. 



6\ = (for n = 1) both satisfy the quartic equation. The 
latter implies that a spin-degenerate zero-energy Landau 
level survives in this context. The two remaining eigen- 
values, e_ and e + , can be estimated perturbatively, in 
the limit A scn huj c <C fi, as e + w — sA so (l — hu c /ti) and 
e_ « —sA so huj c /ti. Their dependence on A so and fku c 
beyond this perturbative regime is shown in the lower 
panel of Fig. 6. This yields the following ordering of 
eigenvalues: = e\ < |e_| < |e+| < |e | = A so . 

The corresponding band structure is shown in the top 
panel of Fig. 6. Its description in terms of spin-polarized 
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bands is slightly more involved than before, but the 
Chern numbers can nevertheless be computed and shown 
to evolve as follows, 



f e T 




= -2 , for - |e | < £ F < -|e+ 


e T 


= o , e 4 = 


-2, for - \e+\ < E F < -\e_\ 


e T 


= i , e+ = 


-2 , for - |e_| < E F < 


e T 


= 2 , e+ = 


-1 , forO<£ F < |e-| 


e T 


= 2 , e+ = 


, for |e_| <E F < \e+\ 


u t 


= 2 , e+ = 


1 , for e + <E F < |e | 



(24) 



indicating that a QH phase is preserved at low energy, 
since the total Chern number never vanishes. This phase 
is peculiar, however, as it is characterized by edge states 
with an unbalanced spin population: |C^| ^ |C^|. For ex- 
ample, two spin-up and a single spin-down counterpropa- 
gating state coexist on the same edge for < Ep < |e_|. 
A given spin species can even become fully gapped, as 
testified by vanishing spin Chern numbers, giving rise to 
spin-polarized edge state transport over a tunable and 
quite large energy window |e+ -e_| ~ A so (l — 2hiu c /ti). 

A QSH phase can be generated close to zero energy by 
lifting the spin-degeneracy of the remaining zero-energy 
Landau level with an arbitrarily small exchange field, de- 
riving from 

H cx = sA cx a r]Q , (25) 

where A ox quantifies the magnitude of the effect. The 
corresponding tight-binding expression is given by 

0i ex = A^J24,a s f Ci,js> ( 26 ) 

i 

using the same notations as before. This yields for the 
spin Chern numbers at low energy: 

C t = -C+ = 1 , for \E F \ < min(|A cx | , | e _| - |A CX |) . 

(27) 

The total Chern number is zero and, contrary to the case 
of layer-symmetric spin-orbit coupling, this time the Z 2 
invariant v = 1, signaling that the QSH phase is non- 
trivial. The energy window where this phase can be 
observed, i.e. the maximum value of the QSH gap, is 
bounded by the value |e_|/2, a lower bound of which is 
given by the perturbative limit Agg^ > A so fkjj c /(2ti), 
as can be checked in the lower panel of Fig. 6. The QSH 
gap could thus potentially reach several tens of meV, al- 
though, in graphene-based systems, it will effectively be 
limited by the highest achievable value of spin splitting 
which should be much smaller [62]. In this respect, and 
despite the need for a perpendicular magnetic field, our 
proposal offers two advantages with respect to that of 
Ref. [35] , where it was recently shown that gated bilayer 
graphene could be turned into a Z 2 topological insulator 
for sufficiently strong Rashba spin-orbit coupling: the 
strength of spin-orbit coupling need not exceed a criti- 
cal value, and spin-orbit coupling need not be present in 
both layers. The latter condition is particularly conve- 
nient if one considers that the most promising chance of 



inducing (intrinsic) spin-orbit coupling in graphene as of 
today is arguably by depositing adatoms on its surface 
[9-11]. 



B. Mechanism II: perpendicular electric field 

Let us now exhibit another mechanism of symmetry- 
breaking which can provide a loophole to circumvent the 
intrinsic difficulty of generating a non-trivial QSH phase 
in bilayer graphene. Forgetting momentarily about spin- 
orbit coupling, let us go back to the Hamiltonian of 
Eq. (14) and consider the effect of an electric field ap- 
plied perpendicularly to the bilayer, 

H v = UaoVz ■ (28) 

This term opens a gap in the energy-momentum disper- 
sion relation by breaking the layer symmetry. It can be 
implemented in a tight-binding model using the following 
expression: 

i£l i£2 

The derivation of the Landau level spectrum requires 
solving once more a quartic equation + a n e^ + f3 n e n + 
7„ = 0, with a non- vanishing linear term /3„ ^ 0: 

a n = - (t 2 + 2U 2 + (2n - 1)A|) 
P n = -2tUA% 

7 „ = U 2 (U 2 + t 2 ) - (2n - 1)A 2 B U 2 + n{n - 1)A% 

Taking into account the spin-degneracy of the spectrum 
and the additional symmetry e„(— r) = —e n (r), one is 
left with two eigenvalues to compute for the lowest en- 
ergy Landau levels, one of which can be easily seen to 
be e+ = tU . The other one, e_, must be computed nu- 
merically. In the limit U,fku c -C t\, it can be estimated 
perturbatively [30] as e_ ss tII(1 — 2hu) c /ti). Its depen- 
dence on U and hio c beyond this perturbative regime is 
shown in the lower panel of Fig. 7. Hence, the QH phase 
has now been gapped by the perpendicular electric field 
at low energy (see top panel of Fig. 7), yielding the fol- 
lowing pattern for the Chern number: 

e t = e+ = -l, for -|e+| <£ F <-|e_| 
C T = 6^ = , for \E F \ < |e_| (31) 
e t = e 4 = l, for |e_| <E F < |e+| 

However, the layer degeneracy of the former zero-energy 
Landau levels has now been lifted, which means that a 
non-trivial QSH phase can once again be generated at 
low energy, provided some spin-degeneracy lifting mecha- 
nism overcomes the gap |e_|. This can be achieved either 
by layer-symmetric spin-orbit coupling (18) or by an ex- 
change term (25). At a critical value of the spin splitting 
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Figure 7: (Color online): (Top panel) Effect of an electric 
field applied perpendicularly to the plane (yielding a layer 
potential asymmetry U = 0.1) on the lowest Landau level of 
bilayer graphene. A gap is opened arising form the lifting of 
the layer degeneracy for the lowest Landau level. Unspecified 
parameter values are the same as in Fig. 5. (Bottom panel): 
Dependence of the eigenvalue e_ on the magnetic field for 
U — 0.01. The dashed line is the analytical prediction (valid 
in the perturbative limit ftw c <C t\) and the thick line is the 
numerical calculation. Inset: Dependence of e_ on the per- 
pendicular electric field for As = 0.05. 



|A CX | = | e | , the lowest bands will cross and give rise to 

a QSH phase 

e t = -e 4 = l, for |^| <min(|A ex |-|e_| , f/-|A ex |), 

characterized by a single pair of counter-propagating 
spin-polarized edge states (see Fig. 8). Once again 
the total Chern number vanishes while the Z2 invari- 
ant v = 1, indicating the non-trivial character of the 
QSH phase. Provided the critical spin splitting could 
be achieved, the maximum value of the QSH gap would 
this time be bounded by the value (U — |e_|)/2, an up- 



Figure 8: (Color online) : Effect of an electric field applied per- 
pendicularly to the plane (yielding a layer potential asymme- 
try U = 0.1) on the lowest Landau level of bilayer graphene, 
in the presence of a spin-splitting term A cx = 0.075. When 
this spin splitting exceeds a critical value, a band crossing 
takes place giving rise to a non-trivial QSH phase. Unspeci- 
fied parameter values are the same as in Fig. 5. 



per bound of which is given by the perturbative limit 
Aqsh < Uhuj c /ti, as can be checked in the lower panel 
of Fig. 7. 

As a closing remark, we note that this second mech- 
anism of symmetry-breaking shares with that demon- 
strated in [35] the property of having edge states in the 
low-energy region of Eq. (32) that are not only spin- 
polarized, but that can also be valley-polarized. This 
can be traced back to the lifting of valley degeneracy 
by the perpendicular electric field, as is apparent in the 
values of the Landau levels given below Eq. (30). This 
valley polarization actually translates into an additional 
topological protection, encoded in the valley Chern in- 
dex v = ^J2 t tG t , with G T = J2 s ^r,s- In the energy 
region of Eq. (32), this index verifies v = 1, indicating 
that the low-energy phase is a so-called quantum val- 
ley Hall phase. The latter is entirely analogous to a QSH 
phase, if one exchanges spin and valley indices: it is char- 
acterized by valley-polarized counter-propagating edge 
states, which can thus only be backscattered by short- 
range (valley-coupling) disorder. Hence, the low-energy 
phase of Eq. (32) should be immune to spin-mixing per- 
turbations as long as valleys remain uncoupled [63] . 



C. Discussion 

Now that we have identified the regimes in which a 
non-trivial QSH phase could arise in bilayer graphene, 
and that we have roughly estimated the order of magni- 
tude of the associated energy gap Aqsh , let us conclude 
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this section by making a few comments on the experimen- 
tal relevance of our results. Until now, we have made the 
natural assumption of disregarding the effect of disorder 
in our system, since one of the essential features of a topo- 
logical phase is its robustness with respect to disorder. 
The presence of the latter could nevertheless prove prob- 
lematic if the typical strength of disorder 5d is ^> Aqsh- 
The available experimental data in graphene-like systems 
seem to indicate that low-energy disorder is dominated by 
charge density fluctuations (electron- hole puddles), but 
the use of BN substrates has been shown to significantly 
reduce their magnitude [36-38]. 

The main other threat to the QSH phase lies in the var- 
ious many-body instabilities which have been predicted 
to occur in bilayer graphene at the Dirac point due to the 
finite density of states. This could lead to a spontaneous 
symmetry breaking of the spin- valley SU(4) symmetry 
in undoped bilayer graphene, causing the emergence of 
a yet unidentified gapped phase, typically of the order 
of a few meV [39-41]. Amusingly, a (many-body driven) 
QSH phase stands among the list of possible candidates 
[42, 43]. 

Estimating the importance of disorder and interaction 
effects eventually boils down to how big a value of the 
QSH gap could be achieved. If Aqsh lies in the 10 meV 
range, then the presence of disorder should be harmless 
to the QSH phase, while actually reducing the effect of 
the Coulomb interaction. On the other hand, if Aqsh is 
rather in the 1 meV range, then chances are great that 
disorder and/or interactions will wash out the picture we 
described. 



V. EXTENSIONS 

Let us now briefly discuss extensions of our model to 
closely related systems. We start by considering different 
types of stacking orders in bilayer graphene, and then 
move on to the case of trilayer graphene. 

A. Other stackings 

The analysis we performed in this article relied on 
the assumption of AB (Bernal) stacking for the bilayer. 
However, other possibilities may occur. One of them is 
the so-called AA-stacking, where both layers are mirror- 
symmetric: A-i atoms sit on top of A\ atoms and B2 
atoms sit on top of B\ atoms. In this case, follow- 
ing the exact same steps as described in section II A, 
the Landau level spectrum can be obtained [44], e^ A = 
±\Ai(ii i |n|fi.w c ). Contrary to the case of Bernal stack- 
ing, the Landau level with lowest energy is no longer 
necessarily that corresponding to n = 0, which leads to 
a peculiar band structure (see Fig. 9) where the low- 
energy physics is described by counter-propagating spin- 
degenerate edge states, characterized by a trivial 6 = 
phase. In a sense, the absence of zero-energy Landau lev- 
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Figure 9: (Color online): Landau level spectrum of 
AA-stacked bilayer graphene with parameter values as in 
Fig. 5. Notice how, at low energy, spin-degenerate counter- 
propagating edge states lead to a trivial topologically trivial 
phase. 



els in this system, which we took as our defining criterium 
for a 2DDFG, is directly responsible for the absence of a 
topological order at zero energy. We additionally checked 
that the symmetry-breaking mechanims investigated in 
this work are ineffective for the present system. 

Besides AB and AA stackings, a whole (continuous) 
family of bilayers referred to as twisted bilayers can be 
studied experimentally. Such bilayers are defined by the 
angle with which the upper layer is twisted from the lower 
layer. This angle can be probed experimentally by char- 
acterizing the induced Moire patterns. Although such 
systems are also interesting in their own right (and ex- 
perimentally relevant), they are not well suited to a tight- 
binding description, especially for small angles, as the 
low-energy physics requires potentially very long-range 
hoppings to be taken into account. We will therefore 
not discuss them any further, and we refer the reader to 
other approaches developed in the literature to address 
their properties (see for example [45]). 

Likewise, so-called double layer graphene [46] - a bi- 
layer where the coupling between the layers is solely 
capacitive (transverse hopping is zero) - crucially re- 
quires electrostatic screening to be taken into account, 
and therefore lies beyond the scope of this paper. 



B. Trilayer graphene 

In the light of our understanding of single layer and bi- 
layer graphene, we finish by briefly discussing how much 
of our previous considerations could find a natural ex- 
tension in trilayer graphene [64] . One generally distin- 
guishes two stacking orders (see Fig. 10): ABA stacking, 
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ABA stacking 



ABC stacking 



A, t D B, A, t 



A: 



Figure 10: Side view of typical stacking sequences of a trilayer 
of graphite: ABA stacking is mirror-symmetric with respect 
to the central layer, while ABC stacking can be seen as the 
natural extension of Bernal stacking in the bilayer (see Fig. 4). 



characterized by a combination of linear and quadratic 
dispersions at low energy, and ABC stacking, character- 
ized by a cubic dispersion and a corresponding diverg- 
ing density of states at low energy which favors many- 
body instabilities. Regardless of the stacking sequence, 
the odd number of layers implies that, as in monolayer 
graphene, including spin-orbit coupling in each layer (via 
a naive extension of Kane and Mele's model) will yield a 
non-trivial QSH phase [65] . The additional presence of a 
perpendicular magnetic field - which has experimentally 
been shown to give rise to a QH effect [47-50] character- 
ized by a spectrum with a 12-fold degenerate zero-energy 
Landau level and described by the following values for 
the Chern number, 

e t = e+ = -3, for -A LL <E F <0 
e t = e l = +3 , forO < E F < A LL l^so-uj 

where All is the energy of the lowest non-zero Landau 
level - will yield a transition from a QH to a non-trivial 
QSH phase at low-energy (top panel of Fig. 11): 



3 , for \E F \ < A s 



C T = G l = 3 , for A so < E F < y/ A LL + A| 



yielding 6 = and v = 1 (mod 2) when \E F \ < A s 



(34) 



Exploring further the fate of the Landau level spec- 
trum, we have checked that applying spin-orbit coupling 
only in the upper layer is (QSH-wise) ineffective. How- 
ever, applying a perpendicular electric field, through the 
tight-binding expression 
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Figure 11: (Color online): Landau level spectrum of (ABC- 
stacked) trilayer graphene with spin-orbit coupling A so = 0.02 
(top panel), and with both perpendicular electric field (U = 
0.1) and spin-splitting A cx = 0.05 (bottom panel). The ef- 
fect of spin-orbit coupling is completely analogous to that 
in bilayer graphene (see Fig. 5), causing a lifting of spin- 
degeneracy in the lowest Landau level. This time, however, 
the odd number of pairs of spin-polarized counter-propagating 
edge states leads to a non-trivial QSH phase at low energy 
(top panel). Additionally, and also as in bilayer graphene, 
the simultaneous presence of a layer-degeneracy lifting elec- 
tric field and a spin-splitting term can also give rise to a non- 
trivial QSH phase at low energy, with a single pair of counter- 
propagating spin-polarized edge states (bottom panel). Once 
more, unspecified parameter values are the same as in Fig. 5. 



%u = -U^c\ci + U^2c\ci, (35) 

has an interesting effect which distinguishes ABA from 
ABC stacking. In the latter case, it opens a gap, while in 
the former it does not (though the QH phase is trivial at 
low energy, due to counter-propagating states) . When an 
exchange term is taken into account, a QSH phase with 
only a single pair of counter-propagating spin-polarized 



states can be accessed (bottom panel of Fig. 11). Thus, 
our second symmetry-breaking mechanism seems to work 
equally well in trilayer graphene, although its relevance 
is debatable in the present context since, as mentioned 
above, a QSH phase could already be obtained in trilayer 
graphene in the absence of any layer inversion symmetry- 
breaking. Additionally, the width of the energy window 
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A 1 = 2 (I 1 lg. 5 top) 


K^ti with c = ±4 
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weak C^brl witn i> = U (mod 2J 


1\ — 2 with A so only in upper layer (.b lg. b top) 


with c = ±1 (spin- unbalanced) 


N = 2 with A so only in upper layer, and A ex 




N = 2 with U (Fig. 7 top) 





N = 2 with U, and |A CX | > |e_| (Fig. 8) 


QSH + QValleyH 


N = 3 


QH with 6 = ±6 


N = 3 with A so (Fig. 11 top) 


QSH 



Table I: Summary of low-energy topological phases in graphene-based 2DDFGs. 



where our mechanism is effective decreases with the num- 
ber of layers, which can be qualitatively understood as 
originating from the proliferation of bands (due to the 
increasing degeneracy of the lowest Landau level) . 

VI. CONCLUSION 

We have considered different examples of graphene- 
based 2DDFGs and shown that, in the presence of both 
spin-orbit coupling and a perpendicular magnetic field, 
a topological phase transition between a QH and a QSH 
phase could take place at low energy. An overall sum- 
mary of the various cases discussed in this article is pro- 
vided in Table I. While the lifting of spin degeneracy in 
the Landau level spectrum was the only requirement to 
observe this transition in monolayer graphene, we showed 
that a similar prescription proves insufficient in bilayer 
graphene, yielding a weak QSH phase at low energy. 

We then proceeded to identify several regimes in which 
a non-trivial QSH phase, characterized by a single pair 
of counter-propagating spin-polarized edge states, can be 
induced in bilayer graphene, all of which involved break- 
ing the layer inversion symmetry. We investigated two 
possible ways of achieving this: (i) by considering the 
presence of spin-orbit coupling only in the upper layer; 



(ii) by applying a perpendicular electric field. In both 
cases, the resulting low-energy phase can then be tuned 
into a non-trivial QSH phase in the presence of an ex- 
change field: in case (i), an arbitrarily small exchange 
term suffices, while in case (ii), a non-zero critical value 
is required. The first of these two cases has the advan- 
tage of requiring only small spin splitting (which will, 
however, effectively control the magnitude of the induced 
QSH gap) and the presence of spin-orbit coupling only 
in the upper layer. The latter condition is crucial, as the 
most promising way to induce sizeable (intrinsic) spin- 
orbit coupling in graphene is arguably by random adatom 
deposition [9-11]. Indeed, the effectively weak intrinsic 
spin-orbit coupling of carbon remains as of today the 
main obstacle in the attempt of experimentally detecting 
the QSH phase in graphene-like systems. In this respect, 
other recently isolated two-dimensional crystals such as 
silicene [51-53] or cold atom optical lattices [54] might 
offer an alternative to probe the physics described in this 
work. 
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